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In this article we give a full description of the dynamics of the flat anisotropic (4+1)- 
' dimensional cosmological model in the presence of both Gauss-Bonnet and Einstein contri- 

butions. This is the first complete description of this model with both terms taken into 
, account. Our data is obtained using the numerical analysis, though, we use analytics to 

explain some features of the results obtained, and the same analytics could be applied to 
higher-dimensional models in higher-order Lovelock corrections. Firstly, we investigate the 
' vacuum model and give a description of all regimes; then, we add a matter source in the 

I , form of a perfect fluid and study the influence the matter exerts upon the dynamics. Thus, 

we give a description of matter regimes as well. Additionally, we demonstrate that the 
presence of matter not only "improves" the situation with a smooth transition between the 
standard singularity and the Kasner regime, but also brings additional regimes and even 

<n : 

Q\ , partially "erases" the boundaries between different regimes inside the same triplet. Finally, 

00 ■ 

^-j- . we discuss the numerical and analytical results obtained and their generalization to the 

CO ' higher-order models. 
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INTRODUCTION 



The idea of extra dimensions traces back to the beginning of last century, to the papers by Nord- 
strom pi, Kaluza ^ and Klein Qj. At that time such ideas were considered as mere mathematical 
speculations, but with time they have come to hold a firm place in the minds of scientists. The 70s 
gave rise to an interest in the extradimcnsional theories due to the development of superstring and 
supergravity theories (see, e.g., [4]). Later in the 90s the interest has only increased owing to the 
possibility to solve the hierarchy problem [5], lower the grand unification scale in the M-theory 6] 



and in many others (see, e.g.,|7|] for recent reviews). 

Dealing with extra dimensions in the cosmological context, one usually uses a modified gravity 
(see the forementioned reviews). One of modifications of this kind is the Lovelock gravity [8|. 
This theory is equivalent to Einstein's theory in (3+1), but starting from (4+1) it gives rise to 
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higher-order (in powers of curvature) corrections. The first of these corrections is known as the 
Gauss-Bonnet term, first found by Lanczos [9[] (therefore it is sometimes referred to as the Lanczos 
term). In the cosmological context the Lovelock (and Gauss-Bonnet in particular) gravity was 
intensively studied over the past two decades 1(3 ]. 

In this paper we investigate the dynamics of the flat anisotropic Universe in Lovelock gravity. 
Since we work in a (4+l)-dimensional case, Lovelock contribution reduces to the Gauss-Bonnet 
term. We took both the General Relativity (GR) and Gauss-Bonnet (GB) contributions into 
account for a reason. Namely, with only one of them taken into account, some exact solutions 



(see, e.g., jlll4l5j]) could be obtained (and that is true not just for GB, but for any order of 
Lovelock corrections and in any number of dimensions [3]). Though when we take into account 
different contributions, the equations become much more complicated making it almost impossible 
to obtain any exact solution in general case (for the first time solutions with both terms taken into 
account were found in [l?]]; in 13] a cosmological solution in exponential form was obtained). Yet, 
the influence of the other (than leading) terms could be significant - for instance, in [ijj], while 
investigating a (4+l)-dimensional model with matter in the form of a perfect fluid in a "pure" 
Gauss-Bonnet case, we discovered an unusual behavior in w < 1/3 case, and we believe it was 
caused by neglecting the GR term in the equations of motion. 



The (4+l)-dimensional case with both GR and GB contributions was studied in 18j], but a full 



analysis of all regimes was not performed there. In this paper we are going to do this - we will give 
a description of all regimes, as well as try to explain some of them analytically. 

The second aim of this paper is to investigate the regimes in presence of matter in the form 
of a perfect fluid. We reported some features of the influence of matter in [l]], but no actual 
description was given. Namely, we reported that the presence of matter leads to an increase of 
the "probability" of smooth transition between the low- and high-energy Kasner regimes. In the 
current paper we give a full description of all types of transitions as well as describe the influence of 
the matter in large. Finally, we generalize the results obtained - both for the vacuum and matter 
cases - to higher even-dimensional models with all possible Lovelock corrections taken into account. 

16]. Indeed, in [lfj we considered the flat anisotropic 



This article is the logical continuation of 



cosmological models in the Lovelock gravity, derived the equations of motion and investigated the 
dynamics with only the highest Lovelock correction taken into account. In this paper we investigate 
the lowest model (D = 4, n = 2) with all corrections taken into account and based on the results 
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predict the behavior of even-dimensional models with all possible Lovelock corrections considered. 



II. EQUATIONS OF MOTION 

Equations of motion can be easily derived from the general ones, obtained in |IG]. hi terms of 
Hubble functions they take a form: dynamical equations 



2(H b + H 2 ) + 2(H C + H 2 ) + 2(H d + H 2 ) + 2H b H c + 2H b H d + 2H c H d + 



+8a 



+ P = 



(H b + H 2 )H c H d + (H c + H 2 )H b H d + (H d + H 2 d )H b H c 

(this is the first equation; the rest of them could be obtained via a cyclic indices transmutation) 
and a constraint equation 



2H a H b + 2H a H c + 2H a H d + 2H b H c + 2H b H d + 2H c H d + 24aH a H b H c H d = p. (2) 

In the vacuum case ([1]) and ([2]) reduce with w = and p = 0. For simplicity we will refer to the 
model with a non-zero matter contribution as the "matter case" (or the "matter model"). 

So in the matter case in addition to ([T]) and ([2]) one also needs an equation of state p = wp for 
the perfect fluid and a continuity equation: 



p+{p + p)(H a + H b + H c + H d ) = 0. (3) 

The system (Hj)— ([2]) in the vacuum and (pQ)-([3]) in the matter case combined with the initial 
conditions completely determines the evolution of the Universe as a whole. Hence we are going 
to solve numerically these equations to find out the future and past evolution of the model with 
particular initial conditions. Scanning over the initial conditions gives us the needed distribution 
of regimes. 



III. VACUUM MODEL: SPECIAL CASES 



Before starting to present the results it would be an asset to consider two special cases that can be 
found on our transition maps. We call the first of them a 3-equal case with three Hubble parameters 



equal to each other. In [18] it was noted that only this regime in a (4+l)-dimensional case has a 
smooth transition between the low- and high-energy Kasner regimes. Using H a = H^ = H c = H 
and H<i = h the dynamical equations reduce to 

6H 2 + AHh + AH + 16aHhH + WaH 3 h + +2h + 2h 2 + 8aH 2 h 2 + 8aH 2 h = 0, (4) 

6H + 12H 2 + 2AaH 2 H + 2AaH 4 = 0, (5) 
and the constraint equation gives h = —H/{1 + AaH 2 ). Solving @ and © gives 



• 2H 2 (l + 2aH 2 ) 

H = l + 4a^ ' (6) 



• 2H 2 {8a 2 H i + 2aH 2 -1) 

(1 + 8a# 2 + 16a 2 //4)(l + 4aF 2 ) ' ^ ^ 



18l |. and we claim ours to be correct 1 . Still, 



Note that our equations are different from those in 
the reason why a smooth transition between the low- and high-energy Kasner regimes occurs is the 
same - the denominator of neither H nor h never crosses zero. 

We call the other case a 2-equal - now only two Hubble parameters are equal to each other. 
Unlike the 3-equal case, here the denominator of Hi can cross zero so a nonstandard singularity 
could occur. But what is important, in this case there also exists a smooth transition between the 
low- and high-energy Kasner regimes, thus, the conditions for a smooth transition between the low- 
and high-energy Kasner regimes are somehow weakened. But since a 2-equal case still needs exact 
equality, the measure of this transition is not improved. 

IV. VACUUM MODEL: RESULTS 

As we claimed above, we are going to make a scans over the initial values of Hubble parameters 



to produce maps of trajectories. Instead of producing 3D maps, as it was done in 



18 ], we will 



1 One can easily verify it - indeed, substituting H a = Hi, — H c = H into Eq. (7) in pj| (that is the dynamical 
equation that corresponds to Hd, so it does not have Hd itself but only H a , Hb, and H c ), one would directly obtain 
([5J; it is quite easy to solve it to obtain ([6]). The difference might be caused by a wrong sign in h = — Hj (1 + 4aH 2 ) 
(in pl| the authors have a plus instead of a minus). 



5 



make several 2D plots with different values for the third Hubble parameter. Thus the 3-equal 
case would be a point on our 2D plot with coordinates that coincide with the fixed third Hubble 
parameter, and the 2-equal case would be a diagonal line and two perpendicular lines with the 
x- and y-coordinates equal to the fixed third Hubble parameter. So we have three "free" Hubble 
parameters, the fourth is calculated from the constraint equation ([2]): 

Tji°) rr(°) _l «"(°) it(0) , rr(0) rr(0) 

Now let us have a look on the denominator of ©. Later, when dealing with the matter case, 
we will see that this denominator (more specifically, its sign) is an essential thing in describing 
the matter regimes. For now let us note that if all three Hubble parameters are positive, the 
denominator is also positive, but if one of the three Hubble parameters is negative, the denominator 
could be (but not always) negative. We refer to the (H^\ H^°\ H^) triplet as a "positive" one 
if it produces the positive denominator of ([8]) and as "negative" if it is negative. 

Thus we fix the initial value for one of the Hubble parameters (without loss of generality, let 



it be H^) and scan over hP' and in [0. . . 1.54- 2] range. In previous studies 12|, ll3j, ll9] we 



usually considered [0, 1] range, but now, for the demonstrative reasons we decided to expand it. 
As for H { a\ we vary it roughly in [—1.5 . . . 1.5] range. Also, as in all previous studies, we consider 
only the initially-expanding Universe: YlnH^ 0. First, we present our results for the positive 
then for the negative. 

Before giving the results, let us summarize what we would expect. For the past evolution, we 
expect either a standard singularity (high-energy Kasner regime) or a nonstandard one. For the 
future, there are three possibilities: a low-energy Kasner regime, a recollapse or a nonstandard 
singularity. For simplicity we denoted them as summarized in Table. With the sign of triplet taken 
into account, we denote the trajectories like their types with appropriate signs, e.g., type VI + 
corresponds to type VI with a positive triplet. 

The results for the > case are given in Fig. [1] In there, we plotted the resulting past and 
future behavior for the fixed Ha , and with vs. as coordinates. The resulting transitions 
are denoted according to the Table, the dashed lines are type I trajectories. The value for Ha 
is increasing from Fig. QJa) to Fig. [Iff): H { a ] = 0.1 in (a), H { a ] = 0.3 in (b), i^ 0) = 0.5 in (c), 
H { a ] = 0.9 in (d), H [ a ] = 0.966 in (e) and H { a ] = 1.1 in (f). 

From Fig. Q]one can understand how the transitions evolve with changing At some value - 



I. 



Table. The classification of possible trajectories in the vacuum (4+l)-dimensional GR+GB model 



From 


To 


Design. 


Standard singularity 


Kasner 


I 




Recollapse 


II 




Non-standard sing. 


III 


Non-standard, sing. 


Kasner 


IV 




Recollapse 


V 




Non-standard sing. 


VI 



in our case with a = 1 it occured at H^ ~ 0.966 - two type-IV regions "detach" from each other; 
with further growth of H a the bottom-left region shrinks; the upper-right one also decreases in 
size and moves by diagonal towards the growth of H^ and Ht ■ The H a ^ value this "detachment" 
occurs at is dependent on a and governed by a high-order equation, so it is impossible to give it 
an exact analytical expression. 

The H { a ] < case is more complicated than the previous one. Indeed, with all three "free" initial 
Hubble parameters being positive, the last Hubble parameter is always negative while ^ H- > 
(one can simply verify this from ([5])). But with being negative 2 the situation becomes more 
complicated. First of all, the calculated Hubble parameter could be either positive or negative. 
First case is partially epistemologically unimportant to us - this case is equivalent to H a > 
and thus it does not bring us new regimes, but for the completeness purposes we keep it. The 
second case could bring new regimes - indeed, with one of the three "free" Hubble parameters 
being negative, the denominator of (jHJ) could be negative and its regimes, as we will see in the case 
with a matter source, are distinct from the positive triplets case. However, two negative initial 
Hubble parameters could also lead to ^ Hf' < 0. 

The results are presented in Fig. [2 Therein we gave 2D scans for four different values of 



H (o). H (0) 



-0.15 in (a) panel, H a 



(0) 



-0.25 in (b), H, 



(0) 



-0.3 in (c) and H a 



(0) 



-0.5 



in (d). At H a ^ < —0.5 type- VI trajectories completely "disappear" from the scene, hence we 
decided to skip the plots of that kind. A white region corresponds to the initial conditions with 
H^ < 0; hyperbola-like regions originate from zeros of the denominator in the Hj expression: 
H { a 0) + H^ ] + H ( c ] + 12<xffl 0) 'H^ 0) 'H ( c ] = 0. Thus, the upper-right "half" of the 2D plots in Fig. [2] 
has a negative denominator of H^ ([8]) while the others have positive denominators. 



2 The same with the previous subsection, we fix Ha^ and plot 2D figs with and as the x- and y-coordinates. 




FIG. 1: Transitions for H a > in the vacuum case. The H a value is increasing from (a) to (f). Trajectory 
types are denoted according to the Table. Dashed lines correspond to the type-I trajectories. 
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FIG. 2: Transitions for H a ' < in the vacuum case. The H a value is decreasing (\H a ' \ is increasing) from 
(a) to (d). Trajectory types are denoted according to the Table. Solid lines correspond to the 1+ trajectories, 
dotted - to the I—. 

In addition to the type II, IV, and VI trajectories at low values of \Hi°^\, one can find type-I 
trajectories as well. There as well exists the transition of I— type, that we can call "the rarest" 
trajectory. Indeed, it exists only at < 0, only at 2-equal lines, at low enough {H^l (type- 
IV should not disappear yet - > — l/(4 v / a)), and on the upper "half" of the plot; the I— 
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trajectories are denoted by a dotted line in Fig. [5]^a). The 1+ transitions are denoted as solid lines 
in Fig. [2]^a) and (b). One can see that type-I trajectories exist between the lower boundary of 
the Yh H i° ] > and 

some point that we denoted as diamonds in Fig. EJa) and (b). These points 
correspond to the 3-equal situation: indeed, from Sec. Ill we remember that in the 3-equal regime 
the calculated Hubble parameter is always negative. If we assume to be that parameter, then 
the 3-equal situation occurs at 



Two different 3-situations are infinitely separated at Ha ~^ (0 — 0), become closer with the growth 
of \H a ' | and coincide when H^ = —l/{<iyfa) (with a = 1 it happened at H a ^ = —0.25; the 
situation is presented in Fig. E^b)). 



V. VACUUM MODEL: DISCUSSION 



In the first half of the article we dealt with the vacuum (4+l)-dimensional flat anisotropic model 
in the Einstein-Gauss-Bonnet gravity. One can see that the type-II trajectories are dominating in 
both H^ > and H^ < situations; in other words, the type-II trajectories are dominating all 
other regimes in the model considered. In the H^ > case we have IV and Vl-type trajectories 
more or less abundant, their presence is in general lessened with the growth of H^\ but they never 
vanish. In the H^ < case the situation is more dramatic: with small \H a \ we have I, II, IV, 
and VI type trajectories, but with the growth of \H a \ we first "loose" I and IV, then VI, therefore 
finally at H^ < — 0.5 (with a = 1; with a different value for a this will happen at another value 
for the H^) we have only type-II transitions. 

In Fig. [3] we presented the examples of all the possible transitions. Black curves correspond to 
the expansion rate in terms of the Hubble parameters Yi Hi > an d grey curves represent the value 
for the denominator of the Hi expression. Time is along x-axis and is normalized in a way to set 
a = 1; zero-point t = corresponds to the start of integration. Type I is in (a), II - (b), IV - 
(c), and VI - (d). Hence there are no type-Ill and V transitions, and now we will try to explain 
why. One can see that the nonstandard singularity occurs only in the vicinity of the standard one 
(but it is not necessary that it will happen - see Fig. [3^a) for example). From Figs. [3]^ a) and (c) 
one can learn that the trajectories, starting from the positive values of the denominator end up in 
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(a) 



(c) 




FIG. 3: Examples of different type of transitions. Black curve is the expansion rate in terms of the Hubble 
parameters y\. Hj, grey curve is the value for the denominator in the Hi expression. Type I is in (a), II - 
(b), IV - (c), and VI - (d). 



the Kasner regime. We checked that for the scans in Fig. rjj and found that the type-IV regime 
completely coincides with an initially positive denominator. The reason behind this link between 
the positivity of denominator and the late-time Kasner behavior is a bit aside from the goals of 
this paper and is yet to uncover - we are going to devote a separate paper to the nonstandard 
singularity in flat anisotropic models in the Lovelock gravity. Anyway, type III and V trajectories 
cannot originate from the positive initial value of the denominator. As for the negative value, the 
remaining two possibilities are presented in Figs. E^b) and (d). One could think of the possibility 
for the type-Ill trajectory to occur on the singular analog of the left branch in Fig. EJb) - on the 
short "window" between the standard and nonstandard singularities, but our numerical analysis 
demonstrated that it is not happening. We believe this situation could not even be constructed: 
crossing the nonstandard singularity changes the sign of the denominator leaving the numerator 
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unchanged; thus at this point the sign of Hi changes, making it useless to try to construct any 
continuation beyond the nonstandard singularity. 



VI. MODEL WITH MATTER 



In this section we start our investigation of the models filled with matter in the form of a perfect 
fluid. In this case the constraint equation gives us a bit different expression for the 4th Hubble 
parameter: 

(0) _ 2fljy + ggggg + 2H^HP - Po 

Prom (jlpp one can easily see that the expression for is split into two parts - p-dependent 
and ^-independent. The latter remains the same with varying p while the former changes the value 
for . At this point the positive and negative triplets start to play an important role. These two 
have a huge difference - the positive triplets have no boundary on the density while the negative 
do. Indeed, keeping in mind that ^ ^ and the fact that the increasing density decreases 
fljj , even if at p = we have Yl^i°^ > 0, at some value for p we will have Yl-^-i = 0- This 
difference between the positive and negative triplets affects the details of their regimes, so while 
describing the regimes we will describe them separately. 



VII. MODEL WITH MATTER: RESULTS 



Now let us describe the matter regimes. The first one to describe is 1+ case . In the vacuum case, 
only 1+ and I— have the smooth transition from the standard singularity to the Kasner expansion, 
but, as we noted in [131 ] . in the model with matter it is no longer the case. In the presence of 
matter the transition remains unchanged: it is from the standard singularity to the expansion. The 
difference between the w < 1/3 and w > 1/3 cases lies in the isotropisation: at w > 1/3 the initial 
singularity becomes isotropic (and the expansion is the GR-dominated Kasner) while at w < 1/3 
we have a standard GB-dominated singularity and an isotropic expansion. 

Trajectories of I— type, as we already mentioned, are "the rarest" trajectories. In the presence 
of matter the regimes changes according to Fig. Ufa); the trajectory types are denoted in the figure. 
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FIG. 4: An example trajectory map in (p,w) coordinates for I— (a), 11+ (b) and II—: w s t < (c) and 
1/3 > w s t > (d) subcases. The stroked region in (a), (c), and (d) panels corresponds to the periodical 
trajectories reported in [13j. See the text for details. 

As a regime with the negative triplet, it has the maximal value for density; the x-axis of the Fig. 
Ufa) covers all possible values for density. 

The transitions in the (p,w) map for 11+ type trajectory are given in Fig. [U(b). There is a 
difference between the regimes at w < 1/3 and w > 1/3: in the w < 1/3 case, the GR-Kasner 
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"transforms" into an isotropic expansion while in the w > 1/3 case, the isotropic one is the initial 
singularity. 

From Fig. H](b) one can easily see what brings an addition of matter to this regime. First of 
all, new, unknown for the vacuum case, type-Ill regime is introduced. It doesn't cover a wide 
range of p and w, but its presence is still important. Secondly, the type-I trajectory in this case is 
unbounded from high densities. Thus, we find that the presence of matter in the form of a perfect 
fluid makes the type-I transition common in terms of initial parameters in contrast with the vacuum 
case, where one requires fine-tuning to achieve this regime. One can also notice that the type-I 
transition covers both w > 1/3 and w < 1/3 cases, so the expansion can be both Kasner-type 
and isotropic, respectively. Finally, the region of type-IV trajectories is also unbounded from high 
densities and lies mostly in the lower half of the (p,w) map; later we will see that this is more or 
less typical for the type-IV trajectories. 

The II— case is more complicated than the others. Namely, it has a different behavior in the 
w s t < and 1/3 > w s t > cases. We consider them separately, and now we only note one feature 
of the w > 1/3 regimes: the initial singularity is GR-dominated in contrast with 1+ and 11+ where 
in the w > 1/3 case the initial singularity was isotropic. 

In w s t < case almost all the trajectories remain the same - II type. Since this is "standard 
singularity —> recollapse" , one can find its "lifetime" . To demonstrate the influence of matter on the 
dynamics, we plotted the resulting lifetimes on the (p,w) map (see Fig. Etc)). The stroked region 
in the bottom-right corner marks the positions of periodic trajectories - a new type of solution in 
the vicinity of the stationary one (p = p s t, w = w s t), described in [ljj. Also one can notice that 
there is no abrupt difference between the w < 1/3 and w > 1/3 regimes. 

And one last point regarding this regime - in the case of w s t < — 1 (this is possible in some 
triplets) we do not have anymore periodical trajectories (and so only type II trajectories remains) 
as they are all located at w < — 1 which is beyond physical assumptions. 

The case with 1/3 > w s t > is different from the described above w s t < - it has more different 
regimes. In Fig. H^d) we presented the map of transitions in (p, w) coordinates. 

The map of trajectories for the IV+ case is given in Fig. [5j There in panel (a) we present a 
large-scale structure of the transitions; in panel (b) - a detailed structure of the boundary between 
I and IV at low densities and in panel (c) - the boundary between I and IV at large. 

The resulting (p,w) transition map for the IV— case is given in Fig. (6]^a). The type IV trajec- 
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FIG. 5: An example trajectory map in the (p,w) coordinates for the IV+ case. In (a) we can see the large- 
scale structure of the transitions (up to very large values of p) , in (b) we presented the boundary between I 
and IV at low densities with feature, and finally in (c) - the boundary between I and IV at large. 



tories are quite bounded in there, but we have another regime, unseen in the vacuum case - type 
V. With it we have all the six possible regimes in the presence of matter. With high enough w s t 
we have a relatively large area with periodic trajectories. 

Two final regimes, VI+ and VI— , are presented in Figs. E^b), (c), and (d). In Fig. [6fb) we 
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FIG. 6: (a) Example of the (p, w) transition map for the IV— case; (b) - the same but for the VI— case. 
One can see the region with periodical trajectories in the bottom-right corner of panels (a) and (b). In (c) 
and (d) panels - (p, w) transition map for the VI+ case: structure at low p in (c) and large-scale structure 
in (d). 



presented type VI— and in Figs. 0(c) and (d) - type VI+: the behavior at low values for p in (c) 
and the large-scale structure of the transitions in (d). 
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FIG. 7: Type-III (a) and type-V (b) transitions in the case with matter. Black line represent JZ-ffi, grey 
denominator of H expression. 

VIII. MODEL WITH MATTER: DISCUSSION 



In the second part of the paper we studied what influence the matter in the form of a perfect 
fluid exerts upon the dynamics of flat (4+l)-dimensional anisotropic models in the Einstein-Gauss- 
Bonnet gravity. One of the two main changes from the vacuum case is a substantial increase of 
the type-I transitions we previously reported in [13]. An interesting feature - only positive triads 
have type-I transition in the presence of matter, negative (with an exception of I—) do not. In this 
paper we demonstrate this increase "quantitatively" . Another feature is the appearance of type-Ill 
and type-V trajectories. First of them appeared exactly the way we discussed in the vacuum case: 
on the singular analog of the left branch of Fig. [3][b). We presented both III and V in Fig. [7p. 
Thus, the presence of matter changes the dynamics to make such trajectories to exist. 

One can find different "features" on the (p, w) transition maps and most of them (those that are 
unrelated to w s t) are just triplet-specific features. A feature at the 11+ transition map (Fig. Hfb)) 
where all five regimes "touch" is a good example. If we draw a w = const line from that point to 
p = we will notice that the lifetimes for type-II trajectories along this line are the same and equal 
to that of the vacuum case. It is obvious that the lifetime in the vacuum case is triplet-specific and 
varies from triplet to triplet, hence the location of features of that type also varies from triplet to 
triplet. The same situation is with the features of other kinds. 



3 Remarks to Fig. ^regarding the time are valid to Fig. [7] as well. 
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Now we want to draw your attention to the II— 1/3 > w s t case. One can see that it is very 
different from other cases. It has periodical trajectories at an "unusual" (from the other regimes 
that have periodical trajectories) place and it has much more type- VI transitions than the other 
regimes do. The reasons behind this behavior are not clear; probably, that is due to crossing 

12 ] numerically for the GB+matter case (later in [3] we 



the "special solution" w = 0, found in 



demonstrated that this is a general feature for all even-dimensional flat anisotropic models in the 
Lovelock gravity with only the higher-order correction taken into account). 

One can notice that generally, inside a triplet, the transition maps have a lot in common. For 
the positive triplet we have III— >I— TV "regimes stratification" (over w) on large p: in the 11+ 
case (Fig. [U^b)), we have the same situation in IV+ (Fig. [5]) and VI+ (Fig. [6]^c, d)). In the 
11+ case type-Ill trajectory "disappears" at some density, but that is just an effect of a particular 
triplet: in another triplet with a lower value for w at which the feature occurs we will have all the 
three regimes at high densities. Negative triplets also have some features in common. They all 
have type-II at the intermediate and high values of density and, if w s t > — 1, they have periodical 
trajectories. They all as well (except II—) have type- VI at low-intermediate values of density. 

Therefore we see that the positive triplets have a high-density regime stratification over w while 
the negative triplets have it over the density itself. Keeping in mind the difference between the 
triplets one can explain this: for negative triplets the density is more important - high density 
can cause Y2i < 0, so the value for density determines initial regimes. For positive triplets 
all values of density are acceptable hence the role of density is less important and the equation of 
state determines the regimes. 



IX. HIGHER-DIMENSIONAL CASES 



Finally, as we claimed at the beginning, we are going to generalize the obtained results on the 
higher-dimensional flat anisotropic models in the Lovelock gravity. First of all, in [lij we obtained 
the equations of motion for the flat anisotropic cosmological model in the Lovelock gravity in 
any number of dimensions with only the highest possible correction taken into account. We also 
demonstrated how one can obtain the equations for the mixture of corrections. If we take into 
account all possible corrections, the resulting equations for different orders would have a similar 
structure, and so would the solutions. Therefore we expect the Il-type trajectories to dominate 
in all the even-dimensional models. The negative multiplets (the analogs of triplets from (4+1)) 



18 



would as well have a more complicated structure than that in Fig. This is caused by more 
possible zeros of the denominator of the calculated Hubble parameter. In its turn, it will increase 
the effective presence of type-I, IV, and VI trajectories. The gain in type-I is also due to the 
increase of possible (D — /c)-equalities, but still the measure will be incomparable with other types 
of regimes. 

We do not expect any changes in the negative multiplets - they are governed by the density - 
but there might be changes in the positive multiplets. Indeed, they are governed by the equation 
of state, and its influence is dependent on the number of dimensions. The equation of state that 
corresponds to the highest possible contribution (n-order) is w eq = l/(2n— 1). With the growth of 
n we have w eq — > (0 + 0); as one of the consequences one would expect a decrease of the presence of 
the II— type trajectory with < w s t < Weq ~ the corresponding range of w would shrink eventually 
decreasing the presence. As another consequence one would expect a decrease of the type-I presence 
in the positive multiplets (apart from I— we do not have type-I in the negative triplets in (4+1)). 
This is caused by the same "shrinking" w eq — > (0 + 0) while w = is the "special solution" for the 
even-dimensional case [3]. Thus the region < w < 1/3 from (4+1) (where the type-I trajectories 
are located) with the growth of n would shrink resulting in a decrease of presence of the type-I 
transitions. 

Also, the characteristic scale of the Hubble parameters where types I, IV, and VI occur will drop 
with the growth of n. This is due solely to the number of dimensions - indeed, the same amount 
of matter produces lower density in the higher dimensions. Our calculations for D = 6 and D = 8 
confirmed that. Thus, purely due to the growth of n the effective presence of the type-I, IV, and 
VI trajectories is decreasing. 



X. CONCLUSIONS 



In this paper we gave the first complete description of all the regimes in the (4+l)-dimensional 
flat anisotropic cosmological model in the Einstein-Gauss-Bonnet gravity. All through the paper we 
hold the condition on the initial values of Hubble parameters: Yl > 0. With the Yl < 
case taken into account, our description will be complete. And the ^ < case could be easily 
obtained from the ^ > case: indeed, assuming ^ < we effectively only make a t — > — t 
transform, so we will have the same, but time-reversed results, such as, "standard singularity — > 
Kasner" turns into "Kasner —> recollapse" and so on. With this remark in mind, our description 
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of the regimes is complete for all the possible values of the initial Hubble parameters. 

We have found that the type-II transition ("standard singularity — > recollapse") dominates 
over all the other types of trajectories in the vacuum case, and describe it both qualitatively 
and quantitatively. We as well found that the type-I transition ("standard singularity — > Kasner 
regime" ) is spread more widely than it was assumed , but still not enough to compete even with 
type IV. 

If we add matter in the form of a perfect fluid, the situation changes drastically. Not only the 
abundance of the type-I transition becomes comparable with abundances of other types [13], but 
we also obtain new regimes that are unknown for the vacuum case 4 (type-Ill and V). With these 
two, in the matter case we have all the possible transitions between two possibilities for the past 
(standard and nonstandard singularities) and three for the future (the Kasner regime, recollapse 
and nonstandard singularity) evolutions. We gave the description of the influence of matter on all 
the initially possible vacuum regimes. 

Apart from the description of the matter regimes, we found that at intermediate-high densities 
all the regimes (inside their triplet) look alike: positive triplets have III— >I— TV "regime stratifica- 
tion" (with a decrease of w) while negative ones have VI— >II "regime sequence" (with an increase 
of p). Thus, the influence of matter is stronger than we anticipated: it does not just significantly 
increase the abundance of the type-I transition and adds the two remaining possible trajectories, 
but also almost "erases the bounds" between different trajectories inside the triplet. But there are 
exceptions like II— with < w s t < 1/3 that we discussed in the relevant section. 

Finally we generalized the obtained results on the higher-dimensional case. We expect that 
the vacuum even-dimensional flat anisotropic models in the Lovelock gravity of the n order are 
dominated by the type-II transition. The presence of the other types of transitions is in general 
more and more suppressed with the growth of n. The presence of matter acts the same (as described 
for (4+1)) way in the case of negative multiplets and the same but with decreases of the presence 
of type-I trajectories with the growth of n in the case of positive multiplets. Therefore, the even- 
dimensional vacuum models are degenerative from the point of view of the presence of a smooth 
transition between the standard singularity and the cosmological expansion; the presence of matter 
in the form of a perfect fluid lifts this degeneration but with the growth of n this improvement 



4 Let us note - considering the ^ < case would not bring these two cases for they are time-reversal to each 
other. 
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decreases. Let us note, though, that (2 + 1) in L\ differs from the described above picture [20]; 
hence, the analysis is valid only at n ^ 2. 
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